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Abstract 

The explicit solutions of the Bogomolny equations for N vortices on a sphere of radius 
> N are not known. In particular, this has prevented the use of the geodesic approximation 
to describe the low energy vortex dynamics. In this paper we introduce an approximate general 
solution of the equations, valid for > N, which has many properties of the true solutions, 
including the same moduli space CP^. Within the framework of the geodesic approximation, 
the metric on the moduli space is then computed to be proportional to the Fubini-Study metric, 
which leads to a complete description of the particle dynamics. 
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1 Introduction 



The abelian Higgs model in the plane is one of the most studied examples of a field theory 
with topological solitons. The solitons are vortices. At critical coupling there are Bogomolny 
equations, and it is known that there is a 2A^-dimensional manifold of gauge-inequivalent A^- 
vortex solutions ||6|. This is known as the A^- vortex moduli space, and denoted Mn- As a 
manifold, Mn = C^. There is a natural metric on A^^r, arising from the kinetic terms in 
the Lagrangian of the model, and it has been proved by Stuart that, at least for finite 
time intervals, geodesic trajectories on the moduli space give a good approximation to the true 
dynamics of slowly moving vortices. 

It is convenient to introduce the standard complex coordinate z on the plane. The locations 
of the vortices are the N unordered points where the Higgs field vanishes. These points may be 
regarded as the roots of a monic polynomial in z (monic means that the coefficient of z^ is 1), 
and the natural coordinates on moduli space are the A^ complex coefficients of such a polynomial. 
If a particular geodesic motion is known, then the time-dependence of the polynomial is known, 
and hence the time-dependence of the roots can be calculated. 

Now, a general formula for the metric on moduli space has been given by Samols |jl2|, but 
it is not explicit, so only very special geodesies, with a high degree of symmetry, are understood 
in detail for A^ > 2 0, ^. One vortex just moves at constant speed along a straight line. 
The geodesic motion for two vortices has been calculated by Samols numerically. The most 
interesting phenomenon is that, in a head-on collision, two vortices scatter at right angles. 
Recently, Manton and Speight |11| have found an explicit metric for A^ well separated vortices, 
from which the geodesic motion could be computed. 

In this paper we are interested in the opposite limit. It is possible to consider the abelian 
Higgs model with fields defined on any compact surface. We shall only consider the case of 
a 2-sphere with its standard round metric, parametrized by its radius R. There are again 
Bogomolny equations and a 2A^-dimensional moduli space of A^-vortex solutions. As a manifold 
this is CP^. However, there is an important geometrical constraint, discovered by Bradlow [Q]. 
This is that the area of the sphere must be greater than 4ttN for non-trivial solutions of the 
Bogomolny equations to exist. Equivalently, > A^. The metric on moduli space is known to 
collapse to zero size as the Bradlow limit R^ \ N is approached. We shall be interested here 
in the case where R^ is slightly greater than A^. One should think of this as a situation where 
the vortices are densely squeezed together. We shall present an approximate general solution 
of the Bogomolny equations, and using this, calculate the metric on moduli space directly from 
its definition. Again the solutions involve a polynomial, and the natural coordinates are the 
complex coefficients of the polynomial. We shall find that the metric is that of Fubini-Study, 
with an overall scale factor that depends on R^ — N. 

The geodesies on Fubini-Study are quite simple, and hence, in principle, the motion of 
vortices can be calculated straightforwardly. However, this does involve finding the roots of 
polynomials with time-varying coefficients, which is not algebraically trivial for three or more 
vortices. We shall present examples, mainly of two vortex motion. We should also remark that 
Stuart's proof of the validity of the geodesic approximation for vortex motion does not extend 
automatically to this regime of being close to the Bradlow limit, so our results on vortex motion 
remain rather formal at this stage. The particle dynamics we obtain at the end is, nevertheless, 
quite "physical". 



1 



The paper is structured as follows. In Section 2 we introduce the Bogomolny equations 
on S^, which is identified with CP^. As in the equations are defined on complex line 

bundles over this surface. In Sections 3 and 4 we explain our approximation for close to 
N ^ and proceed to compute the metric on the moduli space of the approximate solutions. The 
geodesies of this Fubini-Study metric are then used in Section 5 to give an explicit description 
of the A'^- vortex dynamics. Some examples of motions are presented in Section 6, and finally in 
Section 7 a general result concerning the number of vortex collisions is proved. 

2 The Bogomolny equations 

According to Bradlow's generalization of the classical vortices over Q, when the base manifold 
is the sphere 5"^ = CP^ , the setup for the problem is a complex line bundle it : E ^ S'^ oi degree 
N equipped with a hermitian metric h. The Higgs field (p is now a section of this bundle, and 
the gauge potentials are the local 1-forms of an /i-compatible connection D on the bundle. We 
will take the metric on 5^ to he ga := B? x (standard metric on5^) , so that the volume of 
{S^qr) is AttR^. 

Letting A := {/i-compatible connections on E} and T{E) := {global C°° sections oiE} , 
the Bogomolny equations for {D, p) & Ax T{E) are (P, ^) : 

L>°'V = (1) 
F + ^MI-1){voIr) = (2) 

where D^'^ is the anti-holomorphic part of D, vol/j € r2^(S'^, M) is the volume form of the metric 
gji, and —iF is the globally defined curvature form of D, so that F G r2^(5^,M). 

We remark that the problem does not depend essentially on {E,h), because all complex 
line bundles over of a given degree N are isomorphic. In fact, for another choice {E',h') 
there will always be an isomorphism f : E ^ E' such that f*{h') = h . It is then not difficult 
to check that {D,<p) is a solution of (1) and (2) on {E',h') if and only if {f*D , o 0) is 
a solution on {E,h), where f*D is the pull-back connection on E. Notice in particular that 
4> E T{E') and f'^ o 4, e T{E) have the same zeros on S"^, and hence correspond to the same 
vortex configuration. 

We will now define the particular pair {E, h) which is to be used in the remainder of this 
paper. Let Ui = CP^\{[0,1]}, U2 = CP^\{[1,0]}, where we use the standard homogeneous 
coordinates [-20,-21] for points on CP""^, and let : C/j — > C be the standard complex charts of 
CP"*^ with transition functions Lpi o tp'^^ =1^2° <fi^ : z ^ 1/ z . Define Qij : Ui D Uj U{1) by 

521 o V2^{z) = {z/\z\)^, gi2 = l/g2i , gn = 522 = 1 ■ 

Since the gij satisfy the usual cocycle conditions, it is possible to construct a complex line bundle 
TT : E ^ CP^ with triviahzations ipi : TT~^{Ui) ^ Ui x <C such that ipi o ipj'^{p, y) = {p , gij{p) y)- 
The hermitian metric h on E \s defined by requiring the unitarity of the triviahzations Tpi , that 
is \'ip~ {p,y)\\ = |yp ; it is well defined because gij has values in U{1). 
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We should now check that degE = N. Define the real valued 1-forms Ai E ri^(C/j,M) by 

-iN 
2{l + \~z 



Ai = ip*A with A := , , (zdz-zdz) £n^{C,R). (3) 



On Ui n U2 one has 



{^^'r{A,-A2) = (^^Ja -a = I (^A^l" dl^^^"" = i{^^'r{gudg2i) , 
or equivalently 

{-iAi) - (-i^2) = 512^521 , 

which shows that the local forms —iAi and —iA2 define a connection on E. The curvature 
—iF]\f of Dn is a global 2-form on CP"*^ determined by F]y = dAj on Uj. In particular, one can 
compute that 

{vT^TFn = dA = .^-^^^-^dzAdz = (ip^^Yi^vol^) , 

and hence 

degi? := ^ [ i-iFN) = ^ [ {v^'TFn = N . 

Integrating equation (2) over CP"^, and using that f^^i F = 27rdegi? = 27rN , it is clear 
that for R'^ < N the Bogomolny equations have no solution, and that for B?' = N any solution 
[D^cp) must satisfy i;^ = and F = ^vol ^ = F]\f. Since we have already constructed a 
connection Dj^j on E with curvature —iF^, we have an explicit solution of the Bogomolny 
equations for the case = N (which is called the Bradlow limit), and it can be shown that it 
is unique up to gauge transformations. 

For R? > N Bradlow has shown in a more general context, that for any solution (j)) of 
(1) and (2), the section (j) has exactly N zeros (which are called vortices), counting multiplicities. 
Moreover, the moduli space TWat of these solutions up to gauge tranformations, is parametrized 
by the positions in CP^ of these N vortices. Since the vortices are indistinguishable, this moduli 
space is identified with (CP^)^/@Ar, where ©tv is the group of permutations of N elements. 

Now consider the map T : (CP^)^/@7V ^ CP^ defined in homo geneous coordinates by 



0<fc<7V 



[ni,^;!],... ^ [... , ^ . . . . . . 

With some care, one can verify that T is a bijection. In fact, its inverse may be described in 
the following way . Given [wq, ... ,W]\j\ E CP^, consider the non-zero polynomial 

which has degree I < N. Calling zi, . . . ,zi the complex roots of P{z), one has 

T-^[[wo,... ,wn]) = [0,1],... ,[0,1] 

Using this bijection, Mn can also be identified with CP^. 
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3 Vortices near the Bradlow limit 



Although we have an accurate description of the moduh space Ai n , the exphcit form of 
the solutions {D, (p) of (1) and (2) is not known. In particular, this has prevented any 
successful attempt to describe the dynamics of the vortices by means of the well-known geodesic 
approximation. The purpose of this paper is to show that by replacing the exact Bogomolny 
equations by two other conditions, which should be a good approximation near the Bradlow 
limit \ N, one can obtain the solutions explicitly; they also have CP^ as their moduli 
space, and furthermore the dynamics of these "pseudo-vortices" is completely computable in 
the framework of the geodesic approximation. 

Since for B? = N the pair {Dn,0) E Ax T(E) is an exact solution of (1) and (2), we 
may expect that for R"^ close to N the solutions {D, 0) will have D Dn (after a gauge 
transformation if necessary). We therefore impose D = and look for ^ G T{E) such that : 

• (-Dat, (p) is a solution of (1), i.e. D^^cf) = ; (4) 

• {Dn,(P) satisfies (2) "on average", i.e. / ( Fn + - 1) voIr ) = ; (5) 

( We note in passing that eq. (4) is analogous to the equation for electron wavefunctions of the 
first Landau level in the uniform background magnetic field F]\f ; eq.(5) is then a wavefunction 
normalization condition. ) 

We will now find the explicit solutions of (4) and (5), and then describe their moduli space. 
Using the local trivialization ipi of E and the chart ipi of CP^ defined before, the equation 
D^j^<f) = over the domain Ui is the same as {B — iA^'^)(f)i = , or explicitly, using (3): 

dcj>i -Nz ^ 

-d¥ = 2(1 + 

where (pi G C°°(C) is the representative of (j) with respect to ipi . 

Equation (6) has the general solution (f>i = f{z){l + l^p)^^/^ , with / holomorphic on C. 
The section (p oi E determined by (pi, which is only defined over Ui , has a representative (p2 
with respect to ip2 given by (p2iz) = gi2{z)4'i{^) , which is smooth on C\{0}. But since we are 
looking for global solutions of (4), (p must be extensible to all of CP^, and this will happen iff 
(p2iz) is smoothly extensible to C. Writing / as a Taylor series, it is then not difficult to check 
that this requires that the coefficient of z'" vanishes for all n > N. Thus any solution <p of (4) 
must have a representative (pi over Ui oi the form 

aoz^ H \-aN 

and conversely any <pi of this form determines a global section (p of E which is a solution of (4) 
over Ui, and by continuity over all of CP^. 

If (p is represented by <pi as in (7), then the representative (p2 will be 

, / X aoH hoATZ-^ 

= (1 + |,|2)iV/2 
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and hence, as for (1) and (2), any solution cj) of (4) has exactly A'' zeros over CP^, counting 
multiplicities. 

We now turn to condition (5). Using that —iFn is the curvature form of a degree N bundle, 
(5) is equivalent to 

where the last integral is calculated in the appendix for (f>i of the form (7). We can therefore 
conclude that represents a solution (f) of (4) and (5) iff is of the form (7) and satisfies the 
normalization condition 

^ k\{N-ky. , ,2 _ T _ 

Calling D c A X T{E) the subspace of solutions of (4) and (5), we thus get a bijection a : P — »• 
^2A'+i (2 £,^+'^ that maps each cj) ^ V, represented by a 0i like in (7), to the point 

The following step is to determine when two solutions in T> are gauge equivalent. Let 
therefore {Di\f,(j)) and {Dn,(P) be a pair of solutions, an suppose g : CP^ — > U{1) is a gauge 
transformation on E that takes one into the other. Using the usual transformation rule for 
connection forms under g, and the key fact that the connection is fixed, it is readily shown that 
g must be constant. So = e*^(/) for some /3 G M. Since the converse is clear, we conclude that 
{Dn,(P) and (-Dat, (p) are gauge equivalent iff 

^ = e'f^(f> ^ (fi = e^^(^i ^ c = e'f^c 

for some f3 eR, where c, c G 5'2Ar+i ^^.^ ^^le images of {Dj^,(j)) and {DN,(f)) under the bijection 
a. Furthermore, notice that this last condition is also equivalent to 7r(c) = 7r(c) in CP''^, where 
TT : 52^+1 ^ CP^ is the usual principal ?7(l)-bundle. 

Calling Mn the moduli space of solutions of (4) and (5) up to gauge transformations, and 
p : V ^ Mn the natural projection, we therefore have : 



p 



(9) 



where a, defined by the commutativity of the diagram, is, like a, a bijection. The right-hand 
side of this diagram is a concrete model for the space of solutions T> and its moduli space. 
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4 The metric on the moduU space 



Using the usual prescriptions of the geodesic approximation (first described in |jT^ ) , we wih now 
obtain the metric m on A^tv which, within the framework of this approximation, determines 
the dynamics of the "pseudo- vortices" (which from now on wih be just called vortices). 



Suppose one has a curve 7 in P: 

t ^ {DM,m)eV ^ (woit),... ,WN{t))eS'''+K 
A natural hermitian metric on D is defined by 

■■= \ j^^^ K^it), Ht)) voIr , (10) 

where the dot stands for the time derivative. Notice that in this opposed to what 

happens in IQ, the gauge potentials do not contribute to the metric, since the connection in 
our space V is fixed, and thus time independent. 
Writing 

ao(i) -\ h a^it) 

Mt) - (i + |,|2);V/2 

for the usual representative of (j){t), using the unitarity of ^1, and noting that 
^1 = (1 + |2;p)-^/2(do2;^ H h dAr) , one has that 

,(i7 dj^ 1 /■,:,2 2ii22 ^kl{N -ky. . . 

^dt dt'^^' 2 7c (1 + 1-^1) ^ (^ + 1)' 

N 
fc=0 

again using the integral calculated in the appendix. We conclude that the hermitian metric 
on V corresponds via the map a to the restriction to 5'^^''"^ of the canonical hermitian metric 
on C^"^^, up to the constant factor 2'it{B?' — N). This metric will also be called (•,•). 



According to the usual procedure, the metric m on A^jy is induced from (•,•) on T) in 
the following way. Given q ^ V and a tangent vector ^ G TqV , let (^)± be its component 
perpendicular to the subspace of TgD formed by the vectors tangent to curves on D which are 
pure gauge transformations, that is perpendicular to kei{p^,)q. Then 

. * . fdj dj dj d-f 

We will now compute the metric on CP^ corresponding to m by the identification a. It will 
also be called m. 

Using the diagram (9), the subspace keic{p^)q C TgV corresponds to the subspace 
ker(7r=„)„(g) C Ta(q)S'^'^^^ . Given w G S^^^^ C C^"*"^ , we have that ker(7r=K)tt, is the one- 
dimensional real subspace generated by the vector ^(e**i(;)(0) = iw. Therefore, given a tangent 
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vector ^ = (70, . . . ,7^) G T^S^^+^ C T^C^+\ we have 

d-f _ d7 

dt dt [w, w) 

so 

\ ' ' j,k=o 

where the computation in the last step uses that {w,w) = 27r(i?2 - N), since w e 5^^+^ Thus 
TT*m is the restriction to S'^^'^^ of the 2-tensor in 

N 

2tt{R^ — N) {6jk — WjWk) dwj (8) didk 

j,k=0 

Now consider the Kahler form fi associated to m. It is defined, as usual, using the imaginary 
part of m: fi = -Imm e J72(CP^,M). We have 

. N 

TT*// = 7r*(— Imm) = — Im(7r*m) = 2t^{R^ — N)- {5jk — WjWk) dwj A dObk I52JV+1 

j,k=Q 

= 27r(i?2 - iV)l V f , , - ^, du;, A dw^ |^,^+, 

j,k=0 

= 2tt{R'^ - N)^ddlog{\wo\'^ + --- + \wn\'^) \s^n+i = 27r(i?2 - iV) vtVfs 

where //ps is the Fubini-Study form on CP^, and the last equality is a well-known result p. 
160]. Since vr and {tt^)w are both surjective, 

Tr*n = tt*{2tt{R^ - N)fiFs) implies ^l = 2tt{R^ - N)^lFs ■ 

Therefore m = 27r(i?^ — N)mps , where mps is the Fubini-Study metric on CP^, because a 
hermitian metric is uniquely determined by its Kahler form. 

5 Vortex dynamics 

Having determined the metric m on the moduli space Mn = CP^, we will now proceed to 
explicitly describe its geodesies, which provide an approximate description of the low-energy 
particle dynamics. Note that m oc mps implies that the geodesies of m are exactly the Fubini- 
Study geodesies. These are well-known ||8|, p. 277] but nevertheless we will rederive them here 
again. 

Let TT : C^+^ \{0} — > CP^ be the natural projection and xo : f^o — > one of the 
standard charts of CP^, where Uq = {[w^, ... ,w'^]e CP^ : w° / 0}. Cahing {c\ . . . , c^) the 
coordinate functions of this chart, then by definition of the Fubini-Study metric we have on Uq : 

/ipg = laaiog(l + |c^p H h|c^|^) = ^h^i^dc' Adc'' and 

2 2 

"T-FS = hj^ dc' (g) dc'' , 
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with 



h 



■jk 



Ojk 



l + |c|2 (l + |c|2)2 ■ 

For a general Kahler metric the geodesic equations have the simphfied form |14, p. 4]: 



where h^'^ := {sk entry of [h^::]' 



In our case h'^^ = (1 + |c|2)((5*'^ + c^c^), and so the geodesic equations for 
chart xo are 

2 < c , c > . 



1+ < c,c> 



, mps) in the 
(12) 



where c{t) is a curve in and < •, • > is the canonical hermitian product. 
Now consider curves in 5*2^+^ c C^"*"^ of the form: 



7(t) = sin(a;t) y + cos(a;t) x , 



t G 



(13) 



where x,y e C^+i are orthonormal with respect to the canonical hermitian product of C ^ . 
If x^ 7^ 0, then vr o 7(t) ^ Uq only for a discrete set D of non-zero values of t, and a short 
computation shows that, in M\L', 



c(t) := xo o ^ o 7(i) = c(0) + c(0) 



x^ sin{Lot) 



u){y^ sin(ti;t) + x^ cos{ujt)) 



where 



X 



C'(0) = -77 



and 



fc^o 



;(y x' 



x^y^) 



(x0)2 



(14) 



(15) 



One can verify directly that this c{t) satisfies (12), and therefore vr o ^{t) is a geodesic for 
t in ]R\I?, and by continuity for all t. If x*^ = , a similar computation in one of the other 
standard charts of CP^ would establish that, also in this case, vr o ^{t) is a geodesic. 

On the other hand, it is not difficult to verify that every geodesic of (CP'^,mFs) can be 
written as vro^, where 7 has the form (13). Although one could give a general, chart-independent 
argument for this, for later convenience we will proceed unnaturally. Namely, using (15), one 
may simply check that given any c(0) G and c(0) G rc(o)C^ ^ C^, the geodesic ■Ko-f{t) with 



X = (l + |c(0)|2)-i/2 (1, c(0) 
y = ..-Hl + |c(0)p)-V2 ( <m,c{0)> 

u; = (l + |c(0)|2)-l/2 



<c(0),c(0) > 

' C(0) ., , I ,^,|n C(0) 



c(0)p 



l + |c(0)|2 ' 1 + \C 

I < c(o),c(o) > |2y/2 

1 + |C(0)|2 



(16) 



has initial position and velocity Xo ^(c(0)) and (xq ^)*(c(0)), respectively. This shows that every 
geodesic starting in Uq is of the form vr o 'y{t). Similar calculations in the other standard charts 
would extend the result to all of 
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We now note two general properties of the geodesies n o ^(t). Firstly, using (11) and 
(15), one can compute that the velocity of the geodesic, which is a constant of motion, is |a;|. 
Secondly, notice that all the geodesies of {C¥^ ,inFs) are periodic. It is not difficult to show 
that for w 7^ the period is 7r/|a;|. 

We will now use our knowledge of the geodesies on the moduli space (CP^,m) to give an 
explicit description of the vortex dynamics. 

Recall from Section 3 that the solutions {Df^,(j)) G T> with the vortices (zeros of (p) in 
positions ip'^^ {zi) , . . . ,(Pi^{zn) in CPS are represented by a function cpi of the form (7), where 
we now have 

aoz^ H \-aN oc {z — zi) ■ ■ ■ {z — zn) , 

and therefore 

ak = A {-!)'' Sk{zi,... ,zn) , k = 0,...,N 

where the Sj are the usual elementary symmetric polynomials, and A is a normalization factor 
which is non-zero for > N. Thus such solutions {D]\f,(f)) G P correspond by tt o ct to (see (8) ) 

'(-!) ( (^ + 1) ) Sk{z,,...,z^),... 

\ yi\ -t- ij. / J o<fe<Ar 

and by xo o tt o a to c = (c^, . . . , c^) G C-^, where 



c'^ = (-l)'^(^ J Skizi,... ,zn) . (17) 

Inverting these relations, we can obtain the positions of the N vortices as a function of the 
coordinates c'^ of a given point in the moduli space CP^. In particular, to the geodesies c{t) of 
the form (14) corresponds a motion of the vortices determined by: 



k=l ^ ^ 



where the Zi{t) are the coordinates of the vortices in the chart {(fi ,Ui) of CP^. Thus, since 
we know all the geodesies of (CP^,m), we can determine all the possible A'"-vortex motions by 
finding the roots of polynomials of degree A'^ — either analytically for < 4 or numerically for 
N >4. 



Now suppose we are given initial positions ^^(O) and initial velocities ij(0) for the vortices, 
where we assume that the Zj(0) are all different. Through (17) and its derivative we can get the 
corresponding values c(0), c(0) G C^, then use (16) to determine which geodesic c{t) corresponds 
to this initial data, and finally solve (18) to get the motions Zi{t). This general procedure has 
been used to obtain the various special vortex motions shown below. 

We remark that, because (17) is a local diffeomorphism only in the region where the vortices 
do not coincide, only in this region can we guarantee that the final result Zi(t) has indeed the 
prescribed initial velocities. This is why we take the Zi{0) all different. If the Zi{0) are not 
all different there are some values of ii(0) that do not correspond to any Zi{t) coming from a 
geodesic motion. 
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6 Examples of motions 



Using the method described in the previous Section, we now give a few examples of 2-and 3- 
vortex motions on the sphere . The trajectories are shown in the complex plane through the use 
of the stereographic projection ^px : S'^\{A^} — C. The particular initial positions and velocities 
used in each case are listed in Table 1. 





1(a) 


1(b) 


1(c) 


1(d) 


2(a) 


2(b) 


2(c) 


2(d) 


3 


zi(0) 


1 + i 


1 + i 


1 


-2V\/3 


1/2 


1/2 


1/2 


1/2 


a > 


zi(0) 


-1 - i 




-3 


2i/^/3 


i 


i 


i 


i 


i 


Z2(0) 







-1 


l + i/VS 


2 


2 


2 


2 


—a 


Z2(0) 





1 + i 


1 


-l-i/V3 


0.6i 


3.7i 


Ai 


4.5 i 


—i 


Z3(0) 








-l + i/Vs 












Z3(0) 








1 - i/V3 













Table 1: Initial positions and velocities 



To help with the interpretation of the figures, we recall that the stereographic projection 
has the property of mapping circles of S'^ (not necessarily great circles) to circles and straight 
lines in the plane. The inverse (p^^ has the converse property. Also the circle of unit radius is 
always shown in a dashed line; the northern (southern) hemisphere of projects to the exterior 
(interior) of that circle. 




Figure 1 : 

(a) Two colliding vortices, one of which is at rest. 

One of the particles describes a great circle on the sphere that passes through the static 
position of the other. 
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(b) Head-on collision of two vortices with the same speed. 

There are two collisions at antipodal points. The total trajectory is the union of two great 
circles. 

(c) Head-on collision of two vortices with different initial speeds. 

Again two collisions occur. The total trajectory is the union of a great circle and another 
circle. 

(d) Symmetrical collision of three vortices with equal speeds. 

The three vortices collide twice at antipodal points. The total trajectory is the union of 
three great circles. 




Figure 2 : 

Except for 2(c), no collisions occur, and each vortex returns to its initial position after one 
period. The degenerate case 2(c) is the same as 1(c) — one great circle and another circle — in 
a different orientation. 




Figure 3 {with a— 2) : 

No collision takes place and the vortices exchange positions after one period. The coordinate 
c^{t) of expressions (14) and (18) is always zero. The coordinate c^(t) is of the form 
-a^ + 2a(l + a'^)/[2a^ + iuj cot(u;t)] . Taking one of the roots z{t) = x(t) + iy{t) of (18) and 
eliminating t from the system x{t),y{t) , one obtains that the trajectory on the plane has the 
simple equation (x^ + y^)^ + (1/a^ — a^)(x^ — y^) — 1 = 0. These are special cases of Cassini's 
ovals and, when projected back to the sphere, look like the joint of a tennis ball. 
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7 Coincident particles and collisions 



In this Section we start by finding an algebraic condition which determines when a point in 
the moduli space M.j^ = CP'^ corresponds to a vortex configuration where at least two of the 
vortices coincide. We then use this condition to show that, for a system of N vortices starting 
in different positions with arbitrary initial velocities, there are at most 2N — 2 collisions during 
one period of the motion. 

Using diagram (9) and the definition (8) of the bijection a, it is not difficult to recognize that 
a point [wq, ... ,wn] € CP^ corresponds by a.~^ to the equivalence class in A^jv of a solution (j) 
represented by 

where A is a non-zero normalization factor. Thus, asserting that [wq^ ... , wj^] corresponds to 
a solution with at least two coincident vortices is equivalent to saying that one of the following 
conditions holds : 

N 

• P(z) := > TTTTT^ z^~^ has a root of multiplicity at least two ; 

^(/c!(iV-/c)!)V2 

• wq = w\ = , which corresponds to a double zero of (p at [0, 1] G CP^ . 
We now use the following result, whose proof is at the end of this section : 

Lemma: For any n G N, there is a homogeneous polynomial S in n + 1 variables of 
degree 2n — 2, such that S{ao, . . . , an) = if and only if ^^^=0 O'k^^'^ has a multiple root or 
oo = oi = . 

An explicit formula for S is given in the proof and we note that, up to a sign, S coin- 
cides with the discriminant of akZ^"^ whenever ag ^ . 

Using this lemma, it is clear that the points [wq, ... ,wi\[] G CP'^ which correspond to at 
least two coincident vortices, are exactly those of the algebraic hypersurface H of degree 2N — 2 
in CP^ determined by the equation 

As we have seen in Section 5, any motion of N vortices in S'^ corresponds to a Fubini-Study 
geodesic in CP^, and these are all of the form t i— > 7r(sin(a;t)y -|- cos{u}t)x), with x,y G C'^"'"^ 
orthonormal and tt being the projection from C^'^-'^ to CP^. By the discussion above, it is also 
clear that this motion will have a collision of two or more vortices iff the corresponding geodesic 
intersects H. But since this geodesic lies on the projective line L = 7r(spanc{a;, y}) C CP^ , 
and does not intersect itself during one period, we conclude that the number of collisions is not 
bigger than the cardinality of LH H. 

It is, however, a standard result in algebraic geometry that either L G H or # (L n H) < 
degH = 2N — 2 . In fact, denoting x = {xq, . . . ,xn) and y = {yo, . . . ,i/jv) in C^"*"^, it follows 
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that a point tt{ux + vy) in L, with {u,v) G C^\{0}, belongs to H iff 

Q{u,v) := S{uxo + vyo, . . . ,UXN + vi/n) = 0. 

But since S is homogeneous of degree 2N — 2, so is Q, and therefore there is a factorization (see 
i P- 31] ) 

2N-2 

Q{u,v) = {aiU + Piv) , for some (ai,/?^) G . 

If Q is identically zero we have L G H. If Q is not identically zero then the roots of Q are 
{Pi, —ai) 7^ Vi , and it is apparent that Lf] H consists of the points 'K{f3iX — aiy), which are 
at most 2N -2. 

We finally conclude that, either the vortices have a motion with at least two of them 
coincident for all time, which corresponds to L C H, or there are at most 2N — 2 collisions in 
one period. 



Proof of the lemma : 

This lemma is a slight generalization of well-known algebraic results, as stated for example in 
I, p. 168] or I, p. 178]. Consider P{z) = ELo^fc-^""'' and its derivative P'{z) = ELi("' " 
k + 1)0^-12;"'"'^, and form the usual resultant 



TZp^pi{ao, . . . , a„) := 



nao 



nao 



On-l 



where there are n — 1 rows with the coefficients of P and n rows with the coefficients of P' , so 
that the matrix is (2n — 1) x (2n — 1). Applying the usual expansion to the first column of this 
determinant we get 

T^p,P'{ao, . . . , a„) = ao S{ao, . . . , an) 

with 

S{ao, ...,an) = 



ao 



(n — l)ai 
naQ 



ao 

«n-l 



nao 



ai 
ao 



nao 



ao 



nao 
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Expanding again the first columns of the determinants in S{0, oi, . . . , On) we get 

5(0,ai,... ,a„) = (-l)"(n-l)al7^Q,Q/ + (-l)"nal7^Q,Q. = (-l)"(2n - 1) ai 7^Q,Q, 

where Q{z) = Ylk=i '^k z^''^. 

Now, if ao / , by standard results in basic algebra (see 

7^p,p, = (-ir("-i)/2^oi?(P) , 

where D{P) is the discriminant of P, and therefore 

5(ao, • • • , On) = 44> D{P) = <^ P has a multiple root . 
If ao = then P{z) = Q{z) and 

S'(ao,... , On) = 5(0,01,... ,a„) = <^ ai = or TZq^q' = . 

But when oi / 0, by the same algebraic results, TIq,q' = <^ Q = P has a multiple root. 
We finally conclude that 

5(00, . . . , an) = <^ P has a multiple root or oq = oi = . 
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Appendix 



In this appendix we compute integrals of the form 



Ic 



(pip 



(l + |z|2)2 



dz A dz 



(pip 



K2 ^ ^ (1 + a;2 + j/2)2 



dx dy 



where = (1 + \z\^y^/^{aoz^ + • • • + ajv) and V = (1 + |^n-^/2(6o2^ + • • • + 6iv) • 
Write 



N 



(t>tp = ^ aN-jbN-kfkj{z) , 

k,j=0 



with fkj{z) = z^z^l + 



Using polar coordinates and integration by parts, 



I{k,j,N):= [ 



(l + |z|2)2 

27r Sjk 



JO Jo 



oo ^2k+l 



(l+r2)A'+2 



-vr 



(l + r2)W+2 



oo 7 



1 



iV + l7o V (1 + ^2)^+1 



where the vanishing of the boundary terms in the integration by parts is valid for N > k > 1. 
Since 



J(0, 0, iV - A;) = 27r 
we have 

I{Kj,N) = Sjk 



■ dr = 



-TT 



(l + r2)^-'=+2 N-k + l [(l + r2)^-'=+i 



1 



TT 



N-k + l 



k{k-l)---l 



{N + 1)---{N + 2-k) 
The final result is therefore 

?2 N 



mO,N-k) = ^^^^^nS,k 



N 



2iR^ 9 . ^ - ^y^kl(N-k)\ - 

^^^|^|2p dzAdz = 4i?2 2^ l{k, k, N) aN-k hn-k = ^t^R 2^ ^T^-rrf «fc ^fc • 



fe=0 



fe=0 



{N + l)\ 
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